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Abstract 

We derive an explicit matrix representation for the Hamiltonian of 
the Ising model in mutually orthogonal external magnetic helds, us¬ 
ing as basis the eigenstates of a system of non-interacting spin 1/2 
particles in external magnetic helds. We subsequently apply our re¬ 
sults to obtain an analytical expression for the ground state energy 
per spin, to the fourth order in the exchange integral, for the Ising 
model in perpendicular external helds. 


*PACS: 75.10.Pq, 75.10.Jm, 64.60.Cn 
I Corresponding author: adegoke00@gmail.com 


1 



Contents 


1 Introduction 


a 


2 Quantization of a system of non-interacting spin 1/2 par¬ 
ticles in external magnetic fields 

2.1 Change of basis via the eigenstates of the single particle 

Hamiltonian. 

2.2 General basis states for the matrix representation of one 

dimensional spin 1/2 Hamiltonian systems. 


3 


Quantization of the one dimensional spin 1/2 Ising model 
in external magnetic fields 


4 


Example application: ground state energy of weakly inter¬ 
acting spin 1/2 particles in external magnetic fields 

4.1 First order correction to the energy. 

4.2 Second order correction to the energy. 

4.3 Third order correction to the energy . 

4.4 Fonrth order correction to the energy. 

4.5 Approximate analytical expression for the gronnd state en¬ 

ergy per spin for weakly interacting spin 1/2 particles in 
external magnetic fields. 

4.6 Estimation of varions order parameters for the Ising model 

in mntnally orthogonal external magnetic fields. 

4.6.1 Magnetization. 

4.6.2 Nearest neighbonr spin-spin correlation. 


15 

16 

il 

18 

21 


22 


24 

24 

25 


5 Conclusion 


m 


1 Introduction 

Field-induced effects in low-dimensional qnantnm spin systems have been 
studied for a long time IDE]. Hamiltonian models incorporating external 
magnetic fields are gaining popularity among experimentalists as well as 
theoreticians (see references [3112 0 E]). A longitudinal field is often 
introduced mainly to facilitate the calculation of order parameter and 
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associated susceptibility as can be seen for example in references mum, 
and a transverse field to introduce quantum fluctuations [TOl [TT] . 

Our main objective in this paper is to give an explicit matrix representa¬ 
tion for the Hamiltonian of a system of N spin-1/2 particles on a cyclic 
one dimensional lattice chain, interacting via nearest neighbour exchange, 
in the presence of transverse and longitudinal external magnetic helds. 

The Hamiltonian, H, is 

N N N N 

H = -h. St - h, Y S! - A. E E • (1) 

2 = 1 2=1 2 = 1 2 = 1 

where and hy are the uniform external transverse magnetic helds, is 
the uniform longitudinal held, J is the nearest neighbour exchange inter¬ 
action, Si are the usual spin-1/2 operators and the helds hj., hy and hz are 
measured in units where the splitting factor and Bohr magneton are equal 
to unity. Periodic boundary condition is assumed so that = S'/, and 
so on. The parameters hx, hy, hz and J are all assumed to be non-negative. 

It is convenient to write H = Hp + Hj, where 

N 

2=1 

N N N 

■fc, E -'*» E ■5.’-E ■ 

2=1 2=1 2=1 

Hp describes a system of N non-interacting spin 1/2 particles in mutually 
orthogonal external magnetic helds. 

The model ([T]) has been widely studied for various combinations of the 
parameters hx, hy, hz and J, especially for phase transitions (see mu 
[T2] and the references therein). Our aim is to give an explicit matrix 
representation for H, using the eigenstates of Hp as basis. 


Hi = - 
and 

Hp = - 
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Throughout this paper we will make use of the following identities which 
hold for j, k G {0,1}: 

j = sin2(j7r/2), 1 - j = cos^(j7r/2), 

5jk = ^- j - k + 2jk = cos^ {{j - k)7r/2} , 
j + k- 2jk = sin^ {{j - k)7r/2} , 

{-lySjk = l-j-k = 6jk- 2jk = cos{(j + k)Ti/2} , (2) 

in particular (—1)-^ = 1 — 2j = cosjvr, (—1)-^“^ = 2j — 1, 

(-ly + (-1)'= = 2{-iy5,k, = 25,k - 1 = cos{(j - fc)7r}, 

jSjk = jk. 

2 Quantization of a system of non-interacting 
spin 1/2 particles in external magnetic fields 

A system of N non-interacting spin 1/2 particles in mutually orthogonal 
external magnetic fields hy and is described by the Kronecker sum 
Hamiltonian 


Hp^ © Hp^ ® • • • © Hpj^ 


where, for j, fc G {0,1}, each single particle Hamiltonian Hp., at the ith 
site, has the matrix elements, in unit of h. 




with respect to the eigenstates {|Ao), |Ai)} of the spin 1/2 operator Sy 
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whose elements, in unit of h, are 

{Xj\SI |Afc) = cos^ |(j - fc) II = Aj cos^ |(j - ^) 11 • 

The remaining two spin 1/2 operators Sf and S'/ have matrix elements 
given by 

(Ail^riAfc) = ^sin2|(j-fc)|| 
and 

(A^l |Afc) = ^-^l^sin^ i) • 

Parameters hy and hz are the external magnetic helds and a = hx—ihy. 
Explicitly, 

Hf^ = -hxSf - hySf - hzSt 

hz hx thy 
hx T %hy —hz 

2.1 Change of basis via the eigenstates of the single 
particle Hamiltonian 

Solving the eigenvalue equation Hp. \ej) = Sj \ej), the normalized eigen¬ 
states |e:j), j G {0,1}, are found to be 

ki) = acj |Ao) + bjCj |Ai) , 

with corresponding eigenvalues 

Ej =—h/2 cos jir, (3) 
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where 


h — + hi) ^ , 


a = hx — ihy, bj = —h cos jvr — 
and 


(4) 


cos jvr {h + hz cos 

^ (2/i)^/2 {h + hz cos {2hy/‘^bj 

Note that 

a*a = -bobi = hi + = h‘^ - hi, 

a* + a = 2hx, a* — a = 2ihy , 

bjbk = {h + hz cosjTif cos^ f } “ 

^ ^ / cos^ {{j - fc) 7r/2} _ sin^{(j - fc)7r/2} \ 

2h h + h^cosjvr (/i 2 _/^ 2 )i /2 J’ 

and 

1 1 

i=o i=o 



The diagonalizing matrix P has elements Pjk = acfc cos^ (j7r/2) + jbkCk, 
for j,k G {0,1}. Thns, Hp^ is similar to the diagonal matrix D having 
elements Djk = £j cos^ ((j — /c)7r/2), that is 

Hp^ = PDP\ 


P 


Cq(Z Cid \ 

cobo cibi y ’ 


D = 


6 


^0 0 
0 El 






With respect to the new basis, {|£o), ki)}; cind for j, k G {0,1}, the Pauli 
spin matrices have the representation 


{ej\\ek) = cos jTTcos^ |(j - k) || 


+ [{h cos Jti + hz) a + {h cos kir + hz) a*] 


sinMO' - k)7i/2} 
Ah 




{ej\Sf \ek) = cosjTTcos^ |(j - k) || 


and 


+ [{h cos jTT + hz) a 


(hcos /cTT + hz) a*] 


i sin^ {{j — k) vr/2} 
Ah {h‘^-hlf/^ 


hy . 2 ^ 


{ej\S- \ek) = cos jvrcos \ (j - k) 


2h ^ 2 


2.2 General basis states for the matrix representation 
of one dimensional spin 1/2 Hamiltonian systems 

Since Hp is a Hermitian operator that lives in a 2^ dimensional Hilbert 
space, hi, its eigenstates form a complete orthonormal basis, suitable for 
giving matrix representations for operators living in hL and with the same 
conditions at the boundary. The eigenvalue equation for Hp is 

Hp\Er) = Er\Er) , T = 0, 1, 2, . . . , 2^ - 1 . 

For each r the eigenstate \Er) is a direct product of the eigenstates of Elp^ 
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while the eigenvalue is the sum of the respective eigenvalues that is 

N 

l-Zi/j-) l^ri) ® I^J’2) ® ^ l^fi) 

i=l 


and 


N 

Ej. £'^2 d" ■ ■ ■ "h ^rjv ^ ^ ^Ti 1 

i=l 


where 



i = 1,2,...iV, 


where \_z \, the floor of z, is the smallest integer not greater than Thus 
each state \Er) is uniquely represented by a binary vector r = (ri, r 2 ,..., tn)- 


Thus, any operator A inH has the matrix representation A with elements 
given by 

Ars = {Er\ A \Es) . 

Using (IH]) we get 

Er = h^Ti - ^ = hrrir “ ^ • (5) 

i=l 

Note that = Yl!i=i U counts the number of |£i) states in the direct prod¬ 
uct state \Er). The degeneracy of the state \Er) is therefore g^Efl = ^Cmr- 
Thus only the ground state and the most excited state are non-degenerate. 
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3 Quantization of the one dimensional spin 1/2 
Ising model in external magnetic fields 


Explicit matrix representation 


Since Hp is diagonal in the basis the only task is to find the matrix 

elements of H[ and then add them to those of Hp. We have 


N 


H,„ = (-5,1 H, \E,) = -Jj] (5,1 St SI, |5,) 

i=l 

N 


( 6 ) 


lro'_i 


2=1 


where ^ = {Sr^. | Iss^.) and where we have introdnced an At—dimensional 
vector d whose components are 2^ x 2^ symmetric binary matrices d* de¬ 
fined by 


N 


dirs = n 


(7) 


i=i 


Thns = 1 if either the two vectors r and s are one and the same 
vector, that is r = s, or they differ only at the consecntive and {i + iy^ 
entries, otherwise = 0. 


Note that 


dnsidn+isi+idi^^ — driSiCi^s — drs ; 


( 8 ) 


where we have introdnced another At—dimensional vector c whose com¬ 
ponents are 2^ x 2^ symmetric binary matrices Cj with elements given 
by 


N 


Ci^. = 


i=i 


riSi 


(9) 


Thns = 1 if either the two vectors r and s are one and the same vector, 
r = s, or they differ only at the component, otherwise = 0. 
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Motivated by the definitions in Q, ([8]) and ([9]) we introdnce two more 
dimensional vectors, a and (i, whose components are 2^ x 2^ sym¬ 
metric binary matrices, in terms of which the Cj and di matrices may also 
be expressed. The and f3i matrices are defined through their elements 
by 


(yirs = ^nsi = COS^ {(rj - Sj) 7r/2} , 
Ars ~ AiSiAi + l^i + l 


= = cos^ {{vi - Si) 7r/2} cos^ {(rj+i - Sj+i) 7r/2} . 

It is straightforward to verify the following properties for the a* and jdi 
matrices: 

ctiCtj = OLjOti = 2 dijCXi 2 (1 — 

I3i/3j = (3jl3i = 2^ ‘^dijfdi - 1 - (1 — 5ij) 12 ^^ -|- (1 — } 

and 

~ -|- (1 — 5ij) |2'^ -|- (1 — 2^ } ) 

( 10 ) 

where 



(1 1 


J2N = 

1 1 

1 


u 1 

1/ 


is the 2 ^ X 2 ^ all-ones matrix. The ccj and A niatrices are singular and 
have trace equal to 2^. The eigenvalues of a* are 2^“^ repeated twice 
and 0 repeated 2^ — 2 times while those of A 2^“^ repeated four 
times and 0 repeated 2^ — 4 times. Finally using multinomial expansion 
theorem and (ITOj) . it is readily established that the matrices a = 
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and f3 = satisfy 

^2 = 2^-^a + 2^-^N{N - 1) Jsiv, 

^2 ^ 2^-2(a + /3) + 2^-‘^N{N - 3) 

and 


al3 = 2^-V + 2^-^N{N - 2) Jaiv . 

It is now obvious that 

^irs ^rs “1“ (1 ®irs) ^ars,N—l 
= hrs + (1 — C^irs) ^l3rs,N-2 

= Srs + 5prs,N-2 cos^ (aj^,7r/2) , 

^i-rs ~ ^rs "h (1 O^ira'j^i+lra^^rs,!^ — '^ 

+ (1 — ai+i^^)ai^^Si3^^^N-2 

+ (1 - - ai+nJ5/3^^,Ar-3 


= Srs + (5/3^s,Ar-3 + {^l3ra,N-2 “ 5pra,N-z){oiira + “i+lrj 
+ {Sl3ra,N-3 - ‘^Sl3^^^N-2)Oii^a<^i+lra ' 

From ffTT]) and ffT^ we find 


N 

Crs ^ ^ ^ira F <5/3rs,Af—2 

i=l 


and 


N 

drs ^ ^ *^*rs -\- 26 “1“ <5/3rs,V—3 • 

i=l 


( 11 ) 


( 12 ) 
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Explicitly 


cos^ {ai^^TT/2) if - 2 


C?' „ ^ 


0 

1 


if < iV - 2 
if r = s , 


di = < 

Lrs ' 


cos^ {ai^^'n:/2) cos^ (aj+i^^7r/2) 

sinMK.. -«i+i.j7r/2} 

0 


if = AT - 3 

if/3„ = iV-2 
if < AT - 2 

r = s, 


and 


0 



a ^rs<N -2 
ii (]rs = N-2 
ii r = s 


drs ^ 


1 


if < AT - 3 
ii^rs = N -3 


2 ii^rs = N -2 

N if r = s 

V 

From the definitions of the q and di matrices the following additional 
properties are evident; 

1. c" = 2"'“^Cj, for n e Z+. 
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2. The eigenvalues of c* are 0 and 2, each repeated 2^“^ times while 
those of di are 0, repeated 2^ — 2^~‘^ times, and 4, repeated 2^“^ 
times. 

3. The Cj and di matrices are singular and have trace 2^. 

Returning to dH]) and substituting for the matrix elements Sf. , we find, 
after some algebra, 

TT NJh'^z V • 2 r/ \ lc^^ 

Hirs = - i “ ri+l)T^/2] 


„ , ,h.J (ff - ^ 

(1 ^rs) nW) ^rs T (1 ^rs) .i 7_9 ^rs i 


2h^ 


4/i2 


where, (for r ^ s), 


P„ = f;c„cos{(r,_.+n«)^} 

^=1 


TT 


= COS \ {Vk-l + Tk+l) - 


and 


N 


Qrs ^ ^ ^ira') dj-g ^firs, 


2=1 


where 


N 


N 


N 


k = {rj - Sjf = ^J (1 - ^r,s,) = 5^4 si 

j=i j=i j=i 


N-3 


sm" UTj- Sj) ^ 


Explicitly, 


1 a Prs = N - 2 and Vk-i = 0 = r^+i 
Prs= { 0 a Prs < N - 2 or Tk-l + r^+i = 1 

if Tk-i = 1 = rfc+1 
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and 

r -1 ifPrs = N-3 

Qrs \ 

[ 0 iil3rs<N -3 

Putting the results together we hnally have the matrix elements for the 
Ising interaction Hamiltonian, Hi, to be explicitly given by 


Hr = 


NJ hi . Jhl 


N 


4 


+ 2 / 2,2 


i=l 


TT 


sm < [Vi - Tj+i, 


/I ^ ( TT 

(1 - -V.,iV -2 cos I (rfc_i + Tk+l) - 


where 


- (1 - <5r 


J{h^-hl) 

- 


'Jl3rs,N-3 , 


N 

fc = ^jsin2|(r,-s,)^|. 
i=i 


Since H^s = + -^/rs ''^c therefore have that the matrix elements of the 

Ising model in mutually orthogonal external magnetic fields are given by 


rr ^ NJhl^ Jhl^ A . 

Hrs /iuT’s y ^ Ti Oy.^ -|- 07 2 ^ 


2=1 


sm- t (rj - rj+i) ^ 


2 = 1 


/i r r TT 

- (1 - -V.,iV -2 COS I (rfc_i + Tk+l) - 


(1 - 5r-s) 


J{h^-hl) 


4h2 


l3rs,N — 3 5 


with k as defined above. 
Dehning 


/ = h, g = hL_^, F + / = 1, 
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we have 


Hr = 


i=\ 


(1 - ^rs) Vo,iV-2^^ COS I (rfc_i + rfc+i) ^ 


(13) 


- (1 - 6rs)^-Si3rs,N-3 


and 


N - hS 

J^rs TTtrpllOrpg ^ Orpg ^ 0g 


Jf 


N 


^rs ^ ^ ' 


i=l 


IT 


sin < [Vi - Ti+i J 


(14) 


(1 - 5rs)5(3,,,N-2^^^ COS | (rfc_i + Tk+l) ^ 


- (1 - 


,V-3 


where 


= 


N 

E 

i=i 


j’ 


N 

k = ^J 

j=i 


sin 


TT 


Ej - SjJ 


4 Example application: ground state energy 
of weakly interacting spin 1/2 particles in 
external magnetic fields 

When the exchange integral J is small, the Ising interaction term Hi can 
be treated as a pertnrbation oi Hp- In this section, we employ flT^ to hnd 
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corrections, up to the fourth order in J, to the energy of the ground state 
of weakly interacting spin 1/2 particles in mutually orthogonal external 
magnetic fields. Since the ground state of Hp, the unperturbed system, 
is non-degenerate, we will apply the non-degenerate Rayleigh-Schrodinger 
perturbation theory. 


The following particular cases of flTdl) will often be useful. 

p2 7/2 ^ 


Njr Jr 

A ^ 2 




2 = 1 


In particular. 

For s ^ t 


J^ioo - - ■ 


fgJ. 


Hist - - '-^^ 0 st,N -2 cos \ (Sfc_i -f- Sk+l) 2 } ~ ^-^hst,N-3 , 


(15) 


(16) 


(17) 


where 


In particular. 


N 

/c = ^jsin2|(r,-s,)^|. 


i=i 

tilot - - — 0 . 


g^J 


/3ot,N-2 - ^Sl3Qt,N-3 ■ 


Note also from (|5]) that 

Er — Eg = Ers = {mr — 'ms)h, Eqs = —Irish. 


(18) 


(19) 


4.1 First order correction to the energy 

The first order correction to the energy of the ground state of Hp is the 
expectation value of the perturbation Hi in the ground state |i7o) of Hp. 

Thus, quoting ffT6|) . we have 

= {ff/>|E.) = (-Eol H, lEo) = H,„ = -!ffj. (20) 
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4.2 Second order correction to the energy 


The second order correction to the energy of the ground state of Hp is 
given by 


^( 2 ) {Eo\Hj\E:){E,\Hi\Eo) 

° ^ Eo-E, 

s=l * 


|2 

E 


S=1 




According to fflSj) . 


Hi,.= 


fgJ 


Si3os,N-2 


g^J. 


2 4 

42) 


J/3o«,Af-3 ■ 


We therefore see that contributions to Eq come only from states with 
either rris = ^ Sj = 1 (corresponding to Pqs = iV — 2) or ^ s* = 2 

(corresponding to Pqs = iV — 3 in the case when the two |£i) states of the 
direct product state \Es) are consecutive). A typical state with = 1 is 
the state 

\E 2 N- 1 ) = l^i) l^o) ko) • • • ko) • • • ko) = (1) 0, 0, • • • , 0, • • • ,0) 

while a particular state with rris = 2 (and Pqs = iV — 3) is the state 

E 


^ 3 , 2 iv- 2 ) = ki) ki) ko) • • • ko) • • • ko) = (1,1,0, • • •, 0, • • •, 0). 
Therefore 


Hi = and Hi ^ 2 = ■> 

Jq , 2-^-1 2 0 , 3 x 2-^-2 ^ ’ 

and since there are N vectors with = N — 2 and N vectors with 
(3os = N — 3, and using (IT^ . we obtain 


^(2) _ 
U/Q — 


N 


H, 


0,2^^-^ 


N 


Hi 


0,3x2^-^ 


h 


2h 


( 21 ) 


2^2 


Nf^g 


■r 


Ng^ 


r. 


Ah 32h 

The results fl20|) and fl2T]) were also obtained in |13) . 
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4.3 Third order correction to the energy 

The third order correction to the energy of the ground state of Hp is 
obtainable from the formula 




2«_l 2™-l 

EE 

S = 1 t=l 


if, 


h. 




ho 


EoxEi 


H 


ho 




2"-l 

E 

S=1 


\H. 


hs 


^Os 


2^-1 

E 

s=i 


HiJ H, 

El 


2iv_2 2^'^-! 

E E 

s=l £=s+l 


^Iqs ^Ist ^ho 

EosEqi 


H 


ho 


2"-1 

E 

S=1 


if, 


hs 


El 


— Si + S 2 + , 


where 


2^-2 2^-1 


5 *^ = ^ \^hsf S 2 = 2 ^ ^ ^hsHi.tHito 

^Os 


S=1 


5 = 1 £=5 + 1 


Er\f,E( 


Os^Ot 


2"-l 


S 3 = E 


H 


hs 


S = 1 


FS 

^Os 


Note that in the above derivation we made use of the following summation 
identity 

MM M M-1 M 

fst =++ y^ y^ (+++) • 

s=a t=a s=a s=a £=5+1 


Evaluation of Si 

• Contribution from states with = 1 (3qs = N — 2) 

Hios = (from ([HD), + fJ (from ([HD) 

The contribution of the N states with mg = 1 to the sum Si is 
therefore 
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• Contribution from states with = 2 (provided that /Sqs = N — 3 ) 


H 


lOa 



NpJ 

4 


+ PJ 


The N states with = 2, /Sqs = N — 3 therefore contribute 


N 


9^ 

16 



to Si- 


Putting these results together we have 




Npg^J^ 


Ng^J^ 

16 



( 22 ) 


Evaluation of S 2 


S 2 


2^-2 2^-1 


25 : E 

S=1 t=s-\-l 


EosEot 


In each term of the sum, one of four different scenarios is possible, namely, 
rUs = 1 = rrit or rris = 2 = or = 1, = 2 or rris = 2, rrit = 1. We 

look at each possible situation in turn. 


• Contribution to S 2 when rris = 1 = 

In this case, for each s vector, there are two possible t vectors for 
which the matrix element does not vanish, as typihed below: 

s:(0,l,0,0,---,0,0) s: (0,1,0,0,--- ,0,0) 1 ^ _ n 

f : (1,0,0,0,--- ,0,0) f : (0,0,1,0,--- ,0,0) / 

In such a situation. 
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We also have 


|_(s^O,/?o. = iV-2) 


Hi,. = - 
and 

Hi,o = - 7^ 0, Pot = N - 2). 

Since there are N rus = 1 states, the contribution to the sum S 2 
when rris = 1 = is 

(2iV ■ 2 • -fgJ/2 ■ -g^j/4 • -fgj/2)j{-h ■ -2h) 

Npg^J^ 

8h2 

• Contribution to S 2 when = 2 = rrit 

As in the previous case, for each s vector, there are only two possi¬ 
ble t vectors for which the matrix element does not vanish, as 
typihed below: 

s:(l,l,0,0,---,0,0) s: (1,1,0,0,--- ,0,0) 1 ^ . 

f : (1,0,1,0,--- ,0,0) t : (0,1,0,0,--- ,0,1) / 

In such a situation, 

Q^J 

Hist = —^ but hf/.o = hT/o, = 0 since (3to = N - 4. 

There is therefore zero contribution to S 2 when = 2 = rrit. 

• Contribution to S 2 when = 2, = 1 

In this case, typical situations with an s vector and the two t vectors 
for which does not vanish are depicted below 

s:(l,l,0,0,---,0,0) s: (1,1,0,0,--- ,0,0) 1 ^ - 2 

f : (0,1,0,0,--- ,0,0) t : (1,0,0,0,--- ,0,0) / 

From flTT)) we have 

Hi.t = -'^cos(7r/2) = 0, 
signifying a zero contribution to the S 2 sum. 
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• Contribution to S 2 when rUs = 1, = 0 

Here as in the previous case we have = —fgj/2 cos{'k/2) = 0, 
so that again there is zero contribution to the S 2 sum. 

Adding all the contributions we have 


^2 


Npg^J^ 

8h2 


(23) 


Evaluation of 


Q - -u ST 

>J3 -n /00 / . 

s=l 

From (lT6l) . (lT8l) and (IT^ we have immediately that 

53 


(24) 


Finally combining fl2^ . fl23|) and (Elj), we obtain the third order correction 
to the energy of the ground state oi Hp a.s 


( 3 ) mpg^P Npg^P 

° 64/^2 4h2 


(25) 


4.4 Fourth order correction to the energy 


The fourth order correction to the energy of the ground state of Hp is 
given by the standard Rayleigh-Schrodinger perturbation formula 




2iv_i 2 ^^-I 2^^-1 

S=1 t=l U=1 


^ist ^itu ^iuO 

EosEotEou 


H 


loo 


2"_1 2"-l 

EE 

s=l t=l 


HiosHp^Hp^ 


E, 


Os 


El 


- H 


ho 


2«_i 2«_i 

EE 

s=l 4=1 


ElEot 


+ Hl 


2"-l 

E 

S = 1 


H 


hs 


ps 

-^Os 


-E, 


( 2 ) 


2^-1 

E 

S = 1 




hs 


El 


Calculations completely analogous to those in the previous sections, but 
much more involved, give as 

( 4 ) ^ _UNpp 4 55Npg^ . _ Npg^ , _ Ng^ . 

° 192h3 128h3 4h3 2048h3 
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4.5 Approximate analytical expression for the ground 
state energy per spin for weakly interacting spin 1/2 
particles in external magnetic fields 

Adding the energy corrections (j20l) . (El]), (j25|) and (l26l) to the ground state 
energy (obtained by setting = 0 in ([5|)) of the non-interacting spin 1/2 
particles in external magnetic fields we therefore find, to the fourth order 
in the exhange integral, J, that the energy of the ground state, of 

the one dimensional Ising model in mutually orthogonal external magnetic 
fields, for N spin sites is given by 
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where eo = Eq^f/N is the ground state energy per spin, Eq = —h/2 and 
= -J/eq. 

Since p + we can also write 
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or, in a more compact form, 
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Note that when / = 0, then 
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in perfect agreement with the exact result for the ground state energy of 
the transverse held Ising model ra- 
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where is a complete elliptic integral of the second kind. 

The form of fl27|) suggests an exact result for the ground state energy per 
spin of the Ising model in external magnetic helds: 
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where are positive rational numbers, and in particular, Cq™'^ = (— 
for m > 2. 
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4.6 Estimation of various order parameters for the 
Ising model in mutually orthogonal external mag¬ 
netic fields 

The knowledge of eo allows the derivation of approximate analytic expres¬ 
sions for physical quantities such as the magnetization in each direction 
and the spin-spin correlation function for neighbouring spins. 


4.6.1 Magnetization 

Invoking the Hellmann-Feynman rule in ([T]) gives for the x—magnetization 
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TUrr, = 
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dhx dhx dhx h dh^ 


.*=1 / \ Eo ,^) 

and similar expressions for rriy and the y— and z—magnetizations. 
According to (17r|) . 
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SO that for h 7 ^ 0 we obtain 
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Thus for hz < h ^ 0, 
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and for hx < h ^ 0 and hy < h ^ 0, respectively, 
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and 


ruy ^ 


h h ^ r ^ ^ 

I] ^ (2(i + l)f -(m- 1)/) 


m=2 K. k=0 

Note that in the absence of interaction, {z = 0, h ^ 0), + m^ + m^ = 1. 


4.6.2 Nearest neighbour spin-spin correlation 

The spin-spin correlation, is given by 
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yielding 


4 C m—1 

ci,i+i = f+4j2\ E 

m=2 t k=0 

Note that in the absence of interaction, = 0, we have while 

h = hz gives Ci,i+i = 1. 

5 Conclusion 

We have given an explicit matrix representation for the Hamiltonian of 
the Ising model in mutually orthogonal external magnetic fields, with ba¬ 
sis the eigenstates of a system of non-interacting spin 1/2 particles in 
external magnetic fields. We subsequently applied our results to obtain 
an analytical expression for the ground state energy per spin, to the fourth 
order in the exchange integral, for the Ising model in perpendicular exter¬ 
nal helds. Since the Hamiltonian of the non-interacting spin 1/2 particles 
in external magnetic fields is a Hermitian operator that lives in a 2^ 
dimensional Hilbert space, its eigenstates form a complete orthonormal 
basis, suitable for giving matrix representations for any operator living in 
the same Hilbert space and with the same conditions at the boundary. 
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